Gamma Ray Bursts: new rulers to measure the Universe 
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ABSTRACT 

The best measure of the Universe should be done using a standard "ruler" 
at any redshift. Type la Supernovae (SN la) probe the universe up to z~1.5, 
while the Cosmic Microwave Background (CMB) primary anisotropics concern 
basically z ~1000. Apparently, Gamma-Ray Bursts (GRBs) are all but standard 
candles. However, their emission is collimated and the collimation-corrected 
energy correlates tightly with the frequency at which most of the radiation of the 
prompt is emitted, as found by Ghirlanda et al. (2004). Through this correlation 
we can infer the burst energy accurately enough to probe the intermediate redshift 
(z < 10) Universe. Using the best known 15 GRBs we find very encouraging 
results that emphasize the cosmological GRB role. A combined fit with SN 
la yields = 0.37 ± 0.10 and = 0.87 ± 0.23. Assuming in addition a 
flat Universe, the parameters are constrained to be Qm — 0.29 ± 0.04 and fl\ = 
0.71 ±0.05. GRBs accomplish the role of "missing link" between SN la and CMB 
primary anisotropics. They can provide a new insight on the cosmic effects of 
dark energy, complementary to the one supplied by CMB secondary anisotropies 
through the Integrated Sachs Wolfe effect. The unexpected Standard Candle 
cosmological role of GRBs motivates the most optimistic hopes for what can be 
obtained when the GRB-dedicated satellite, Swift, will be launched. 

Subject headings: Gamma Rays: bursts — Cosmology: observations 

1. Introduction 

Recently, Ghirlanda, Ghisellini and Lazzati (2004, GGL04 thereafter) found a surpris- 
ingly tight correlation between the peak of the 7-ray spectrum E peak (in a vF v plot) and 
the collimation corrected energy emitted in 7-rays E 1 in long Gamma- Ray Bursts (GRBs). 
The latter is related to the isotropically equivalent energy -E 7i i SO by the value of the jet 
aperture angle 9, by E 1 = E^ tiso (l — cos 6). The scatter around this correlation is tight 
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enough to wonder if the correlation itself can be used for a reliable estimate of E 1yiso , making 
GRBs distance indicators, and therefore probes for the determination of the cosmological 
Qm, parameters, and for the exploration of the matter to vacuum dominance transition. 

This issue is similar to the case of SN la: they are not perfect standard candles (i.e. 
their luminosities are not all the same), nevertheless the luminosity of a specific supernova 
can be found through the correlation of their luminosity and the speed of the decay of their 
light curve (i.e. the slower the brighter, Phillips 1993; Riess, Press & Kirshner 1995). It is 
the existence of this correlation among SN la which made possible their cosmological use 
(Riess et al. 2004, hereafter R04; Perlmutter et al. 1999; Schmidt et al. 1998). 

Very recently, this problem has been explored by Dai, Liang & Xu (2004), which found 
tight constraints on Qm and using the correlation found by GGL04. Their result is 
however based on a strong assumption: they assume as universal the correlation measured 
in a particular cosmology (without errors on its slope), and use it to derive the cosmology 
itself. Actually, the best fit correlation depends on the cosmology adopted to derive burst 
luminosities and the correlation should be re-calibrated for each cosmology. 

In this letter we demonstrate that a correct approach leads to less tight constraints 
on the cosmological parameters using GRBs alone. On the other hand, a more interesting 
cosmological result can be acheived if a combined fit with SN la is performed. 



As in the case of SN la, the use of a class of objects as cosmological "rulers" requires 
them being standard candles. The luminosity of GRBs, calculated assuming isotropy, spans 
~ 4 orders of magnitude (Frail et al. 2001), but strong observational evidence (i.e. the 
achromatic break in the afterglow light-curve) indicates that the burst emission is collimated 
into a cone/jet of some aperture angle 6 (Levinson & Eichler 1993; Rhoads 1997; Sari, Piran & 
Halpern 1999; Fruchter et al. 1999). The corresponding energy emitted in 7-rays, corrected 
by the collimation factor (1 — cos#), clusters around E 1 ~ 10 51 erg, with a small dispersion 
(0.5 dex), yet not small enough for a cosmological use (Bloom et al. 2003). 

GGL04 found a tight correlation between E 1 and the (rest frame) peak energy E veak 
of the vF v prompt emission spectrum: E 1 oc E* cak . The exact value of x depends on the 
assumed cosmology. Using Vt M = 0.3 and H A = 0.7 we have 

/ 7-, \ 1.416±0.09 



The scatter of the data points around the correlation is of the order of 0.1 dex. This allows 



2. The Hubble diagram of GRBs 
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to reconstruct the value of E 1 by measuring -Z? pea k- 

This is analogous to SN la, for which there is a tight relation between their peak lu- 
minosity and the stretching factor of their optical light curve (Phillips 1993; Goldhaber et 
al. 2001), with less luminous supernovae showing a faster post-maximum light curve decay 
(Reiss et al. 1995). The proper modelling of this effect (Hamuy et al. 1996; Perlmutter et 
al. 1999) improves the determination of the SN la luminosity and consequently reduces the 
scatter in the Hubble diagram, yielding constraints on the cosmological parameters (see R04 
using SN la with redshift up to z ~1.75). 

The E 1 — -Speak correlation for GRBs makes them a new class of "rulers" for observational 
cosmology and combining GRBs and SN la can further reduce the region of allowed values 
in the cosmological parameter space. Furthermore, since GRBs are detectable at larger z, 
they are a powerful tool to explore in more detail the cosmic kinematics. 

The difference between the standard candle assumption and the use of the intrinsic cor- 
relations, for both GRBs and SN la, is shown in Fig. 1 (top and bottom panel, respectively) 
through the Hubble diagram in the form of luminosity-distance vs redshift. In the upper 
panel we assume that GRBs and SN la are standard candles with a unique energy (E 1 = 10 51 
erg) for GRBs and luminosity (B—-21.1) for SN. The derivation of the luminosity distance 
D L for SN follows straightforwardly from their distance modulus (R04). For GRBs we have 
D L = (1 + z)E 1 /[AnJ : 1 (l — cos#)], where JF 7 is the 7-ray fluence (i.e. the time integrated 
7-ray flux). Note that the determination of 9 requires the knowledge of the isotropic energy 
(see e.g. Eq. 1 in Frail et al. 2001), in turn requiring specific values of (^m, ^a)- In the 
bottom panel we plot SN la and GRBs after correcting for the stretching-luminosity and the 
E^-E^k relations, respectively. In this case the isotropic energy E^^ so of GRBs has been 
estimated from their measured -E pea k through the E^-E^k correlation and the error on the 
slope of this correlation has been properly included in the total uncertainty. Also in this 
case we must fix a given (Qm, ^a) cosmology both for the derivation of 9 and for finding 
the best E~ r E pcak relation. As in the SN la case, the luminosity distance of GRBs derived 
from their E 7 -E pcak correlation (bottom panel) highly reduces the scatter around possible 
different cosmologies (solid, dashed and dotted lines). Moreover, GRBs populate the z > 1 
region, where D L (z) is rather sensitive on (Qm, ^a)- 

3. Constraints on cosmological parameters 

The correlation found by GGL04 assumes f2 M =0.3 and il\=0.7, and h = 0.7. Changes 
on Q M and Qa induce a change on the normalization and slope of Eq. 1, together with a 
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different scatter of the data points around the best fit line. We can then ask what is the 
pair of cosmological values Qm, which produces the "minimum scatter" around the fit, 
performed using the very same Qm, pair. To this aim we use all the 15 bursts of known 
redshifts, E peak and jet break time tbrcak listed in Tab. 1 and Tab. 2 of GGL04. 

The difference with the study of Dai, Liang & Xu (2004) lies mainly in this point: they 
assumed that the E^-E peak correlation is exact and cosmology independent, while it is not 1 . 

Additional differences concern: i) the estimate of the errors on the density of the inter- 
stellar medium when it is unknown (they assume n = 3 ± 0.33 cm -3 while we allow n to 
cover the entire 1-10 cm -3 range); ii) we do not exclude GRB 990510 and GRB 030226 from 
the analysis; iii) we include GRB 030429, for which a jet break time was recently found by 
Jakobsson et al. (2004); iv) we always use (1 — cos 9) (instead of the 8 2 /2 approximation) as 
the collimation correction factor (also when estimating the error on E^j. 

We also consider the 156 SN la of the "Gold" sample of R04, finding the corresponding 
Qm, ^a contours using the distance moduli and corresponding errors listed in their Tab. 5. 
Fig. 2 shows our results. GRBs alone (red lines) are almost insensitive to Q\ but limit Qm to 
lie within ~0.05 and 0.22 (68% confidence level). 

We also show the region pinpointed by the WMAP experiment (Spergel et al. 2003), 
which is only marginally consistent with the allowed region from SN la alone (blue lines in 
Fig. 2). The combined GRBs+SN la fit (filled regions in Fig. 2) selects a region which 
is more consistent with the cosmic microwave background (CMB) results. The minimum 
(with a reduced Xred =1 - UQ ) is for ^m=0.37±0.15 and ft A =0.87±0.23 (1-a). Assuming a 
flat universe yields fi M =0.29±0.04 and fi A =0.71±0.05. 

If we use the "classical" Hubble diagram method, we compare the D-^ values given by 
estimating E 1 through the actual correlation found in each point of the Qm, ^a plane with 
the luminosity distance calculated through e.g. Eq. 11 of R04 (see also Carrol et al. 1992). 
Then, by a x 2 statistics, we find the confidence regions in the Qm, ^a plane, which are 
plotted as dashed line on Fig. 2. This classical method is very similar to the previous one, 
since it uses the same available information. The small difference (contours slightly larger) 
is due to the fact that with the "minimum scatter" method we calculate the distance of 
the data points from the correlation (i.e. perpendicular to the fitting line), while, using the 
"classical" Hubble diagram method, we are using the distance between the E 1 data point 



iFor instance, using ft M =0.4, fi A =0.6 results in £ 7 = (3.7 ± 0.9) x 10 50 (£ pcak /267 keV) 1 - 38±0 09 erg 
(i.e. a ~2.6% and ~16% change in slope and normalization with respect to Eq. 1). With Om=1, Ha=1 
we obtain £ 7 = (3.0 ± 0.9) x 10 50 (£ pcak /267 keV) 129±0 - 08 erg (i.e. a -9% and -30% change in slope and 
normalization with respect to Eq. 1). 
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and the corresponding E 1 by the correlation. 

We can further constrain, with the combined GRB and SN samples, the dark energy 
component which is parametrized by its equation of state P = wpc 2 . Furthermore, w could 
be varying, and one possible parametrization is w = wo + w'z (see e.g. R04). Adopting 
this law, we compute the luminosity distance according to Eq. 14 of R04, assuming a flat 
cosmology with Q M =0.27. In this case the fit is performed in the w -w' plane for GRBs, 
SN and for the combined samples. As before, we recompute the E^-E pea k relation for each 
Wq, w' pair 2 . Fig. 3 reports the corresponding confidence intervals. Besides making the 
confidence region smaller than what derived for SN alone, the effect of GRBs is to include 
within the 68% contour of the joint SN+GRB sample (filled region) the wo — — 1, w' — 
point, corresponding to the classical cosmological constant. 



4. Discussion 

GRBs can now be used as cosmological "rulers" , bridging the gap between the relatively 
nearby type la supernovae and the cosmic microwave background. The SWIFT satellite 
(Gehrels et al. 2004), designed for the fast localization of GRBs, is expected to find about 
one hundred of GRBs per year: it can open up a new era for the accurate measurements of 
the geometry and kinematics of our Universe (for a more extended discussion see Firmani 
et al. 2004). We stress that, besides finding high redshift bursts, which are of course very 
important for finding tighter constraints, it is crucial to find low redshift GRBs, to determine 
the i? 7 -i? P eak correlation in a redshift range which is not affected by the Qm, values. This 
would allow to use the resulting correlation unchanged for all values of Qm, ^a, strongly 
reducing the associated errors. In turns, this will allow to constrain cosmological parameters 
independently from SN la. This is important since GRBs are unaffected by dust extinction 
and it is very unlikely that two completely different classes of objects would have similar 
evolutions to mimic a consistent set of cosmological parameters. 

In Fig. 4 we show an illustrative example of what can be done if a given correlation were 
known to be valid independently of the cosmological parameters. For this we have chosen the 
correlation given by Eq. 1. It can be seen that even the limited sample of our bursts can 
strongly influence the GRB+SN confidence contours, making them more in agreement with 
the WMAP results (not unexpectedly, since we have used just the correlation appropriate 



2 As an example of how the correlation is sensitive to the change of w , w', consider that, for w = —0.7 
and to' = 0.2 the correlation becomes E 7 = (3.75 ± 0.90) x 10 50 (£ , poak /267 keV) 1 - 37±0 - 09 erg (i.e. a -3.4% 
and ~15% change in slope and normalization with respect to Eq. 1). 
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for n M =0.3 and fl\=0.7. A similar consideration concerns the Dai et al. 2004 result). 

We would like to stress that in order to use GRBs to find the cosmological parameters, 
we need a set of well measured data, and especially a well measured jet break time tbreak, 
necessary to find the collimation angle 9, and a good spectral determination of the prompt 
emission. 

We thank Annalisa Celotti for useful discussions. DL thanks the Osservatorio Astro- 
nomico di Brera for the kind hospitality during part of the preparation of this work. G. 
Ghirlanda thanks the MIUR for COFIN grant 2003020775D02. 
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Fig. 1. — Classical Hubble diagram in the form of luminosity-distance vs redshift z for 
Supernova la (open green circles: Calan/Tololo sample (Hamuy et al. 1996); open blu circles: 
Perlmutter et al. 1999) and GRBs (filled red circles: the 15 bursts in GGL04). In the top 
panel the SN la and GRBs are treated as standard candles (no corrections applied); for GRBs 
E 1 = 10 51 erg is assumed. In the bottom panel we have applied the stretching-luminosity 
and the E^-Ep^ relations to SN la and GRBs, respectively, as explained in the text. Note 
that, for GRBs, the applied correction depends upon the specific assumed cosmology: here 
for simplicity we assume the standard f2M=0.3, f2A=0.7 cosmology. Both panels also show 
different D L (z) curves, as labelled. 
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Fig. 2. — Constraints in the Om-Oa plane derived for our GRB sample (15 objects, red 
contours); the "Gold" Supernova la sample of R04 (156 objects, blue contours, derived 
assuming a fixed value of H Q = 65 km s _1 Mpc -1 , making the contours slightly different 
from Fig. 8 of R04). The WMAP satellite constraints (black contours, Spergel et al. 2003) 
are also shown. The three colored ellipsoids are the confidence regions (dark green: 68%; 
green: 90%; light green: 99%) for the combined fit of SN la and our GRB sample. For GRBs 
only, the miminum Xred = 1-04, is at f2M=0.07, Q\=1.2. 
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Fig. 3. — Constraints on the w , w' parameters entering the equation of state p = (wq + 
w'z)pc 2 , where p is the dark energy density, wo = — 1 and w' = correspond to the 
cosmological constant We assume a flat geometry and f2M=0.27 (see also R04). Blue 
contours: constraints from type la SN (R04). Red contours: constraints from our GRBs, 
Colored regions: combined constraints (dark green, green and light green for the 68%, 90% 
and 99% confidence levels, respectively). 
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Fig. 4. — Example of how GRBs can contribute to the determination of the cosmological 
parameters once the E 1 — -E pea k correlation will be found in a cosmology independent way 
(i.e. finding bursts at small redshifts). For this example we assume that the correlation 
of Eq. 1 is valid for any cosmological parameter. We show the contours in both the Qm, 
Qa plane (main figure) and in the wo — w' plane (insert, a flat cosmology with Qm=0.27 is 
assumed). Lines and colors are as in Fig. 2 and Fig. 3. 



